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We study the classical dynamics in a generic first-order quantum phase transition between the 
U(5) and SU(3) limits of the interacting boson model. The dynamics is chaotic, of Henon-Heiles 
type, in the spherical phase and is regular, yet sensitive to local degeneracies, in the deformed phase. 
Both types of dynamics persist in the coexistence region resulting in a divided phase space. 
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The interacting boson model (IBM) [l| describes quadrupole collective states in even-even nuclei in terms of a system 
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of N monopole (s) and quadrupole (d) bosons, representing valence nucleon pairs. In addition to its traditional role of 
interpreting spectroscopic data, the model provides a fertile ground for studying quantum phase transitions (QPTs) 
and mixed regular/chaotic dynamics in a mesoscopic (finite) system. QPTs refer to structural changes induced 
by a change of parameters A in the quantum Hamiltonian, H(\). The underlying mean- field (Landau) potential, 
V(X), determines the nature of the QPT. In particular, for discontinuous (first-order) QPTs, V(X) develops multiple 
minima that coexist in a range of A values and cross at the critical point, A = A c . In the IBM, the integrable dynamical 
symmetry (DS) limits relate to stable structural phases and QPTs are obtained by mixing terms from different DS 



chains [2j|. The competing interactions that drive these QPTs can affect dramatically the nature of the dynamics and, 
| in some cases, lead to the emergence of quantum chaos. 

QPTs 0d3| and chaos have been studied extensively in the IBM, albeit, with a simplified Hamiltonian which, 
I . for first-order QPTs, gave rise to extremely low barrier and narrow coexistence region. Recently, the identification of 
IBM Hamiltonians without such restrictions [3] , enabled a comprehensive analysis of generic first-order QPTs between 
spherical and deformed shapes fcllioj]. The dynamics inside the phase coexistence region was found to exhibit a very 
simple pattern. A classical analysis revealed a robustly regular dynamics confined to the deformed region and well 
separated from a chaotic dynamics ascribed to the spherical region. This divided phase space structure manifests 
itself also in the quantum analysis, disclosing regular rotational bands in the deformed region amidst a complicated 
environment. In the present contribution, we illuminate the origin of this intricate interplay of order and chaos for a 
— . particular case, where the stable spherical and deformed phases posses U(5) and SU(3) DS, respectively. 
^ \ Apart from kinetic rotational terms, the relevant Hamiltonian (up to a scale) is 

•i-H ■ 

#i(p) = 2(1 - 2p 2 )h d (n d - 1) + 2R\{p) ■ R 2 {p) , (1) 
H 2 {0 = ZPlPb + Pl'Pi ■ (2) 

Here n d is the d-boson number operator, R^(p) = v^s^t + py/7(Sd^)f\ P f = d) ■ S - 2(s^f , P^ = 2std+ + 
y/7(cft<ft)p and the centered dot implies a scalar product. The parameters that control the QPT are p and £, with 
< p < l/y/2 and £ > 0. For p = and £ = 1, one recovers the U(5) and SU(3) DS limits, where the Hamiltonians 

ffi(p = 0) = 2h d [2N-n d -l] , (3) 
H 2 (Z = 1) = [-C su(3) + 2N(2N + 3)] , (4) 

involve the relevant Casimir operators. The two Hamiltonians of Eqs. (U)-© coincide at the critical point p c = l/V2 
and£c = 0: H 1 (p c )=H 2 {£ >c ). 

The classical limit is obtained through the use of coherent states, rescaling and taking N — > oo, with 1/N playing 
the role of h. The derived classical Hamiltonians, Hi(p), %(£), involve complicated expressions of shape variables 
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(a) p = 0.03 (b) p = 0.2 ( c ) p = .2 (d) p c = 0.707 




FIG. 1: Poincare sections (top row), plotted at y — 0, and selected trajectories (bottom row), depicting the classical dynamics 
of Hx{p) © for several values of p and energies E. (a) p = 0.03, E = Ex/21, (b) p = 0.2, E = 5Ex/12. (c) p = 0.2, 
E = 12_Bi/21. (d) p c = 1/V2, E = £i/12 (the critical point). The energies are given with respect to the domain boundary 
Ex =V x {pr,P = V§) =2. 



(/3, 7), Euler angles and their conjugate momenta. Setting the latter to zero, yields the following classical potentials 



V x (p) = 2f3 2 - /3 4 /2 - 2p ^/? 3 cos 3 7 , (5) 
Va(0 = £ (4 - 6/3 2 + 9/3 4 /4) + 2(3 2 - (3 4 /2 - ^4 - 2/3 2 /3 3 cos3 7 . (6) 

Vi(p) and V2(0 serve as the Landau potentials with the equilibrium deformations (/3 cq , 7 oq ) as order parameters. The 
potential Vx(p) [Va(£)] has a global spherical [deformed] minimum with, respectively, /3 oq = [/3 cq = 2/V3, 7 eq = 0] . 
At the spinodal point (p* = 1/2), Vi(p) develops an additional local deformed minimum. The two minima become 
degenerate at the critical point p c = \j\[2 (or £ c =0), and are separated by a barrier of height Vb — 0.268 (compared 
to Vb — 0.0018 in previous works [3]). The spherical minimum turns local in V^(^) for ^ > ^ c and disappears at 
the anti-spinodal point (£** = 1/3). The order parameter /3 eq , is a double- valued function in the coexistence region 
(in-between p* and £**) and a step-function outside it. 

The classical analysis simplifies considerably when the dynamics is restricted to L = vibrations. In this case, the 
classical Hamiltonians, Hx(p), become two-dimensional in the polar coordinates (3 g [0, y/2], 7 € [0, 27r) and 

momenta pp £ [0, y/2], p 7 G [0, 1] (or, equivalently, in Cartesian coordinates, x = /3cos7, y = f3smj and p x , p y ). The 



classical motion can then be depicted conveniently via Poincare surfaces of section ll| , shown at prescribed energies 
in Figs. [T][21 along with selected trajectories. 

The classical dynamics on the spherical side of the QPT (0 < p < p c ), is governed by Hx(p)- In the U(5) limit 
(p = 0), the system is integrable and Uxifi = 0) = 2(T + /3 2 ) - (T + /3 2 ) 2 /2, where T = pj + /3' 2 p 2 . As shown 
in Fig. HJa), the sections, for small p, show the phase space portrait typical of a weakly perturbed anharmonic 
(quartic) oscillator with two major regular islands and quasi-periodic trajectories. The effect of increasing p on the 
dynamics in the vicinity of the spherical minimum (x « 0) , can be inferred from a small /3-expansion of the potential, 



Vx(p) ~ 2f3 2 — 2pv / 2/3 3 cos 37. To this order, V\{p) coincides with the well-known Henon-Heiles (HH) potential 12 1. 
As shown for p = 0.2, at low energy [Fig. [ljb)], the dynamics remains regular, and two additional islands show up. 
At higher energy [Fig. Hfc)], one observes a marked onset of chaos and an ergodic domain. The energy where the 
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FIG. 2: Same as in Fig. [T]but for the classical dynamics of -ff2(£) © for selected values of £ and energy E — E2/2I, where 
E2 = V-2.{£; ft = V2) = 2+£. Panel (a) shows a typical pattern encountered for most values of £ > 0. Panels (b)-(c)-(d) portray 
the surfaces in the presence of local degeneracies of normal modes, ep/e~, = 1/2, 2/3, 1, respectively. Panel (d) corresponds to 
the SU(3) DS limit, Eq. ©. 



order-to-chaos transition occurs decreases with p. This typical HH-type of behavior persists in the vicinity of the 
spherical minimum throughout the coexistence region, including the critical point [Fig. [ljd)]. It is also present in the 
region where the spherical minimum is only local (0 < ^ < £**), since V^(£) ~ 4£ + (2 — 6£)/3 2 — 2ft 3 cos 37. 

The classical dynamics on the deformed side of the QPT (£ c < £ < 1), is governed by and has a very different 

character, being robustly regular in the vicinity of the deformed minimum (x ~ 1). At low energy, the motion reflects 
the ft and 7 normal mode oscillations about the minimum. As shown in Fig. [2ja), the trajectories form a single 
set of concentric loops around a single stable (elliptic) fixed point. They portray 7-vibrations at the center of the 
surface (p x 0) and ft- vibrations at the perimeter (large \p x \)- This regular pattern of the dynamics is found for most 
values of £ > both inside and outside the phase coexistence region. Noticeable exceptions occur in the presence 
of resonances, which appear when the ratio of normal mode frequencies, R = e^/e 7 , is a rational number. At low 
energy, this happens at discrete values of the control parameter £ ss £r, in a narrow interval around £r = (3i? — l)/2. 
Panels (b)-(c)-(d) of Fig. [2] show examples of such scenario for R = 1/2, 2/3, 1. The corresponding surfaces exhibit 
four, three and two KAM islands, respectively. The phase space portrait for (£ = 1,R = 1), shown in Fig. [2jd), 
corresponds to the integrable SU(3) DS limit, Eq. ((4]), and is the same for any energy. Interestingly, in the coexistence 
region, where the Landau potential accommodates both the spherical and deformed minima, each minimum preserves 
its own characteristic dynamics. This is evident in Fig. [U(d) and Fig. EJa), which depict the dynamics of the same 
system [Hi(p c ) = #2(6=)] at the same energy, but in different regions (i-ranges) of phase space. 

In summary, we have shown the distinct morphology of classical orbits and different onset of regularity and chaos 
in a generic (high-barrier) first-order QPT between the U(5) [spherical] and SU(3) [deformed] phases of the IBM. 
The spherical phase displays a chaos-susceptible dynamics, similar to the Henon-Heiles system, due to a ft 3 cos 37 
perturbation, significant near the spherical minimum. In contrast, the low energy dynamics in the deformed phase is 
robustly ordered, and reflects the ft-j vibrations about the deformed minimum. This regular dynamics is sensitive to 
local degeneracies of these normal modes. The two types of dynamics preserve their identity and can be detected as 
long as the relevant minimum exists. This leads to a divided phase space structure throughout the phase-coexistence 
region, as found in fsl-tlol|. 
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